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We investigate dynamics of (4+1) and (5+1) dimensional fiat anisotropic Universe filled 
by a perfect fluid in the Gauss-Bonnet gravity. An analytical solutions valid for particular 

> 

' values of the equation of state parameter w — have been found. For other values of 

w structure of cosmological singularity have been studied numerically. We found that for 
w > 1/3 the singularity is isotropic. Several important differences between (4+1) and (5+1) 
dimensional cases are discussed. 

I. INTRODUCTION 



X 

. Recent several years show increasing interest to modified gravity both for (3 + l)-dimensional 

Universe (which is mostly motivated by attempts to explain observed accelerated expansion of our 
Universe) and for multidimensional cosmological models which are more speculative, though they 
pose important questions on the nature of cosmological singularity and possible evolution of the 
Universe in its very early epoch. In the latter perspective the Lovelock gravity [ijl is one of the 
ver y p opular extension of General Relativity (GR) (see, for example, 

12, 13, Q, [3]). It keeps the order of corresponding equations of motion unchanged with respect 



to GR, and , as it was recently claimed in nice thermodynamical properties of GR survive in 
Lovelock gravity in contrast to, for example, f{R) theories. 

In the present paper we consider (4 + 1) and (5 + 1) dimensional anisotropic Universe. In 
these dimensions the only non-Einstein term in the Lovelock action is the famous Gauss-Bonnet 
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combination. Power-law solutions which replace Kasner regime in the Gauss-Bonnet gravity in a 
flat anisotropic multidimensional Universe have been studied since the end of 80s [l^]- However, all 
previous works have been devoted to vacuum solutions. Importance of the Kasner solution (which 
is a vacuum anisotropic solution in GR) follows from the fact that an ordinary barotropic fluid 
different from a massive scalar field is dynamically unimportant near a cosmological singularity. It 
means that we can still use the Kasner solution to describe initial singularity in a Universe filled 
by rather general kind of an ordinary matter. The goal of the present paper is to provide similar 
analysis for multidimensional Universe in the Lovelock gravity. 



II. MAIN EQUATIONS 

We consider a flat anisotropic metric in (n-l-l)-dimensional space-time. We are dealing with 
Einstein-Gauss-Bonnet gravity, and non-vacuum space-time. Lagrangian of this theory have a 
form 



C = R + aC2 + Cm J 
where R is Ricci scalar. Cm is the Lagrangian of matter fields and £2 



C2 = R^ua/sR"'"''^ - ^R^.uR'"' + R^ (1) 

is the Gauss-Bonnet Lagrangian. 

In the present paper we are interesting in the behavior mostly in the vicinity of the cosmological 
singularity, which allows us to take into account only corrections of the highest possible order. 



In our case it is the Gauss-Bonnet contribution, so we ne; 
matter sources this problem have been studied in 
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gleet Einstein terms. In the absence of 



18l ]. in the present paper we take matter 



into account. The full Einstein - Gauss-Bonnet system shows a complicated behavior even in the 
vacuum case [l^, [2^ , and we leave investigation of such system with matter for a future work. 
We are working in the flat background, so the metric we considering has a form 



g^,y = diag{-l, aj{t),al{t),...,al{t)}. (2) 

We use perfect fluid with the equation of state p = wp as a matter source; after varying action 
obtained from (HD using the metric above and perfect fluid as a matter field, one can obtain the 
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following equation of motion: there are n dynamical equations and a constraint equation. The _7-th 
dynamical equation has the form 



- 3 Yl HiHkHiHm - Y,{Hf + Hi) ^ HkHi = map; (3) 

i,k,l,mj^j 
i<k<l<m 

while constraint equation becomes 



i<k<l<m k<l 



3 J2 HiHjHkHi = ap, (4) 

i<j<k<l 

where a = 2'!TG/ac^ is a combination of fundamental constants. We assume a > in order to avoid 
instabilities in quantum theory. Also we need the continuity equation: 



p + ip + p)J2Hi = 0. (5) 

i=l 

To demonstrate the difference between (4+1) and (5+1) cases let us write down respective 
equation separately for these two cases. For (4+1) case we have 

{Hb + Hi)HcHd + {He + H^)HbHd + {Ha + Hj)HbHc + wap = (6) 

as a dynamical equation (the rest three can be obtained by cyclic transmutation of indices) and 
the following constraint: 

SHaH^HcHd = ap. (7) 

For (5+1) case we have 

SHhH^HaHf + {Hh + H^){H,Ha + H^Hf + HaHf) + {H^ + Hl){H^Hd + H^Hf + H^Hf) 

(8) 

+{Hd + Hj){HbHe + HbHf + H^Hf) + {Hf + H)){HbH, + Hi,Hd + H,Hd) + = 
as dynamical equations (rest four can be obtained the same way as for (4+1) case) and 
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3 [HaHf)HcHfi + HaHfyHcHf + HaH^HdHf + HaHcH^Hf + HbHcH^Hf] — ap 



(9) 



as the constraint. 

Now one can clearly see the differences between these two cases. First of all, the structure 
of the dynamical equations in the (4+1) case ([6]) lacks the term with the product of four Hubble 
parameters. Indeed, this term in the j-th equation consists of the product of four Hubble parameters 
with indices different from j and each other. In (4+1) dimensions we have only three indices different 
from a given j, so this term is absent. In the (5+1) case ^ there exits exactly four such indices, 
so one possible term exists. In higher dimensions with a bigger number of Hubble parameters 
the corresponding sum has more than one term. The second important difference exists in the 
structure of the constraint equations - in the (4+1) case the left-hand side of ([7|) has only one term 
(all possible sums of four different Hubble parameters needed) while in the (5+1) case ^ it has 
five of them. Both these features result in interesting differences in dynamical behavior of (4+1)- 
and (5+l)-dimensional models (see below). 



Now let us find some exact solutions with one particular anzatz. These solutions are of power- 
law- type, so scale factors take a form ai{t) = 1^% after substituting Hubble functions Hi = Pi/t and 



III. POWER-LAW SOLUTIONS 



their derivatives Hi 



Pi/t^ into ^ and (HD one can rewrite constraint equation 
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(10) 



i<j<k<l 



as well as equations of motion 



3 ^ PiPkPm 



'm 





i,k,l,mjtj 
i<k<l<m 



i^j k,l^j,i 
k<l 



there pQ is the matter density at some given moment of time. 

It is clear that if p 7^ then these equations can be solved only under condition 




(12) 



5 



Hence, field equations are 

3 ^ PiPjPkPl = crpo, (13) 
i<j<k<l 

-3 ^ PiPkPiPm-^Pi{Pi-l) ^PkPi = wapQ. (14) 

i,k,l,m^j ijtj k,l^j,i 

i<k<l<Tn k<l 

In Einstein gravity Kasner solution (for vacuum) implies YliPi ~ such a condition 

conserves in Jacobs solution (for maximally stiff fluid) (for Kasner and Jacobs solutions see e. g. 
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22(1) ■ Assume now that similarly in Gauss-Bonnet gravity condition ^^Pi = 3 (obtained earlier 
for vacuum) conserves for some kind of matter. One should keep in mind, though, that for vacuum 
(4+1) case we still have YliPi ~ 3 but under condition that one of pi is always zero - indeed, using 
anzatz above to Eq. ([7]) one can see that one of pi should be always zero. 

From X^jPi = 3 and (fT2]) one can obtain w = 1/3. Translate now equations (fT3l) and (fHl) while 
taking w = 1/3 and J2iPi — 3 into account. For example, 

24 ^ PiPjPkPi = ^Pi^Pj J2 P'' J2 

i<j<k<l i jy^i ky^ij l^i,j,k 

= '^P^'^Pj Pki^-Pi-Pj -Pk) = 

= ^YjP'Y.P^ UPk-^Ypi^Pj YPk = ■■■ = 

i j^i k^i,j i jy^i ky^ij ^ 

= 81-54j;p? + 24^p3_g^^4+/'^^A . 

i i i \ i / 



(15) 



At that all equations ()13l) . (I14p turns out to be identical with each other. Any of them takes a form 
81-54j]p2^24^p3_g^p4^ /"^pA ^8^^^_ (-^g) 

i i i \ i / 

Therefore, with u; = 1/3 considered metric with pointed above anzatz will be an exact solution 
only if 

Y^Pi = 1^^ =^^ 3 ^ PiPjPkPi = o-po- (17) 

i i<j<k<l 

Cases n = 4 and n = 5 were investigated for arbitrary w. At n = 4 there are no other anisotropic 
power-law solutions with nonzero energy density, but for n = 5 there is the following solution: 

1 

Pl=P2=P3=P, Pa=P5 = -^P, 
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p-2 3 4 

w = , apo = --P . 

p 4 

This solution have a remarkable feature: if 3 visible dimensions expanse isotropically then extra 
dimensions contract isotropically. However, here apo < 0, therefore, since a = either po < 0, 
or a < 0. The latter means no global energy minimum and hence instability in quantum theory. 
So this solution seems to have no physical meaning; we wrote it down here just for mathematical 
completeness. 



IV. NUMERICAL ANALYSIS 



In order to understand what happens for w 7^ 1/3 we remind a reader similar situation in the 
Einstein gravity. Kasner solution gives us a condition ^^Pi = 1 which indicates that the volume 
expands linearly with time, and, correspondingly, the matter density in the vicinity of Kasner 
regime decreases as p ~ y-w-i ^ ^-w-i other hand, left-hand side of equations of motion 

decreases as These two scaling rates coincide at w = 1 (the Jacobs solution). For w < 1 a 
small amount of matter becomes dynamically unimportant for evolution towards a cosmological 
singularity, and ultimately disturbs the Kasner solution for an expanding Universe. 

Of course, in the Einstein gravity we know the situation much better, because the dynamical 
equations describing Bianchi-I Universe with a perfect fluid can be solved exactly. This solution 
tells us that the Kasner regime is the unique past attractor in the case w < 1 independently of 
the initial conditions, while isotropic regime is the unique future attractor [22]. Though w can not 
exceed the value of 1 for any stable "physical" matter, there are some situation (for example in 
brane cosmology [Q]) where the dynamics can be described as it would be driven by a matter with 
an effective w > 1, and in this case the picture is the opposite one with past isotropic and future 
Kasner attractors 
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25|. 



In the Gauss-Bonnet gravity the situation is as follows. Left-hand sides of Eqs. (3, 4) scales as 
in vacuum, while volume increases as t^. Adding a matter we obtain p ~ V^^^^ ~ ^-Stu-s j£ 
the density of matter is small enough to keep the volume expansion rate close to its vacuum value. 
Comparing these two rates we see that the power-law vacuum solution is destroyed by matter if 
w > 1/3 for a contracting Universe and for w < 1/3 during expansion. On the other hand, the 
power-law solution is stable with respect to adding a small amount of matter with w < 1/3 for 
contracting and with w > 1/3 for expanding Universe. 

In order to study global stability of the power-law solution and check whether its instability 
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FIG. 1: Behavior of individual Hubble parameters (a, d, g), Kasner exponents (b, e, h) and relative Hubble 
differences (c, f, i) for [(5+1), w < 1/3] (first row), [(5+1), w > 1/3] (middle row) and [(4+1), w < 1/3] 
(last row) cases. Note that a;-axis t is not a cosmological time, but rather an internal parameter of our 
numerical calculations. This implies for both Fig. 1 and Fig. 2. 

for appropriate values of w leads to isotropisation we provide numerical analysis of the system (3, 
4) for various initial conditions and w. The results we have found are different for (4 + 1)- and 
(5 + l)-dimensions. 
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FIG. 2: Behavior of individual Hubble parameters (a, d), Kasner exponents (b, e) and relative Hubble 
differences (c, f) for [(4+1), w = —0.01] (first row), [(4+1), w > 1/3] (last row) cases. 

In the (5 + l)-dimensional case the situation is similar to the General Relativity dynamics 
described above. In Fig. 1 (a, b, c) we plot typical behavior of power indices for w < 1/3. There we 
have individual Hubble parameters in (a) panel, individual Kasner exponents in (b) and relative 
differences {Hi — Hi where i ^ 1) m (c) panel. We have found that for all studied initial conditions 
the past singularity is anisotropic and described by above mentioned power-law solution. And 
although in Fig. 1 (b) four Kasner exponents tends to the same value (which is equal to zero), 
from Fig. 1 (c) one can verify that the singularity is anisotropic - differences are non-zero and not 
equal to each other. On the other hand, the singularity is isotropic (again, independently of initial 
conditions) if it; > 1/3 (see Fig. 1 (d, e, f)). In this case one can see that all Kasner exponents tend 
to the same and non-zero value (Fig. 1 (e)) and Hubble differences all tend to zero (Fig. 1 (f)). 

The situation in (4 + 1) dimensions is more interesting. First of all, Eq. (6) indicates that while 
in the vacuum case at least one Hubble parameter should vanish, even a small amount of matter 
makes solutions with vanishing Hubble parameters impossible. This means that an arbitrary small 
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p changes properties of the solution significantly and can not be treated as a perturbation at any 
given moment of time. Thus, dynamics of (4 + 1) dimensional Bianchi-I Universe can show features 
completely different from those in the vacuum case. 

We have found such qualitative difference in case of pressureless matter w = 0. It is easy to 
see that the conditions pi = P2 = pi = Pa = ^ satisfy the system (5, 6) (wc use pi = —Hf/Hi), 
providing an exact solution, which describes anisotropic Universe with constant (given by initial 
data and not changing with time) ratio of Hubble parameters. Note that for a higher-dimensional 
Universe a product of four Hubble parameters in the Eq.(5) destroys this solution, so it is valid 
only in (4+1) dimensions. 

Typical behavior of power indices for other w < 1/3 is shown in Fig. 1 (g, h, i). It is quite 
similar to the (5 + 1) case with u; < 1/3, but there are some differences, like tending some Hubble 
parameters to zero. In theories with Einstein term in the action this feature can violate dominance 
of Gauss-Bonnet contribution near a singularity, so this case needs more work. Initial singularity 
appears to be anisotropic with at least one of the Hubble parameters tending to zero. The last 
property is absent for w = 0, and this value is an exceptional one: any nonzero w leads to a 
qualitatively similar pictures, as it can be seen in Fig. 2 (a, b, c) plotted for w = —0.01. Since 
w = —0.01 satisfy w < 1/3 condition, it looks the same with "usual" [(4+1), w < 1/3] case. With 
w ^ the sudden drop of Kasner exponents to zero becomes sharper and sharper, while in limit 
■u; = it disappears, and all Kasner exponents become equal to unity all the time, as explained 
above. 

On the other hand, we have found that independently of initial conditions and the equation of 
state parameter, the initial singularity for the case of u; > 1/3 is isotropic (see Fig. 2 (d, e, f)). The 
situation is the same as for [(5+1), w > 1/3] case. 

V. CONCLUSIONS 

We have studied the cosmological dynamics of a flat anisotropic multidimensional universe filled 
with a barotropic fluid in the Gauss-Bonnet gravity. Power-law analytical solutions existing for 
particular value of the equation of state parameter (w = 1/3) which are analogous to the Jacobs 
solution in the General Relativity existing for w = 1 have been found. 

In the (5 + l)-dimensional case this analogy exists also for other values of w, leading to isotropic 
nature of the initial singularity for > 1/3 and anisotropic singularity for li; < 1/3. In (4 + 1) 
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dimensions the picture is less clear. First of all, an exceptional solution with constant anisotropy 
exists for w = 1. Second, though our numerical results show clearly that singularity for < 1/3 is 
anisotropic, some Hubble parameters tend to zero, indicating that neglected Einstein term in the 
action may be important in this case. This issue requires special investigation. In case of u; > 1/3 
singularity is isotropic in (4 + 1) as well as in (5 + 1) dimensions. 
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